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CN ■ Abstract 

> : 

■ We extend the Kolmogorov phenomenology for the scahng of energy spectra in high-Reynolds 



number turbulence, to explicitly include the effect of helicity. There exists a time-scale th for 
helicity transfer in homogeneous, isotropic turbulence with helicity. We arrive at this timescale 
using the phenomenological arguments used by Kraichnan to derive the timescale te for energy 
^ transfer (J. Fluid Mech. 47, 525-535 (1971)). We show that in general th may not be neglected 

^ \ compared to te, even for rather low relative helicity. We then deduce an inertial range joint cascade 

^ ' of energy and helicity in which the dynamics are dominated by te in the low wavenumbers with both 

energy and helicity spectra scaling as k~^/^; and by th at larger wavenumbers with spectra scaling 
as /c~^/^. We demonstrate how, within this phenomenology, the commonly observed "bottleneck" 
in the energy spectrum might be explained. We derive a wavenumber which is less than the 
Kolmogorov dissipation wavenumber, at which both energy and helicity cascades terminate due to 
dissipation effects. Data from direct numerical simulations are used to check our predictions. 

PACS numbers: 47.27.Gs,47.27.Jv,47.27.Eq 
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Energy and helicity are the two known inviscid invariants of the Navier-Stokes 

equations. It was postulated in that in isotropic flows with helicity, these quantities 
cascade together from large to small scales. This joint forward cascade of energy and helicity 
has been verified by direct numerical simulations, most recently at a resolution of 512^ 
gridpoints ^. Kraichnan defined the shear time-scale te for energy transfer, based solely 
on energy dynamics. Assuming that helicity dynamics are also controlled by te-, a k~^^^ 
inertial range scaling was established for both energy and helicity spectra P]. 

We would first like to ascribe spatial geometrical properties to the types of quantities 
used to derive the relevant timescales. We recall the spectral formulation ('Uj(k)'u*(k)) of 
the two-point velocity correlation function in isotropic, homogeneous, statistically stationary 
turbulence. It may be decomposed into its index-symmetric and index-antisymmetric parts 
as 

E,,(k) = i((S.(k)n*(k)) + (n,(k)S*(k))), (1) 
4(k) = ^((^,(k)«*(k)) - («,(k)ti*(k))), (2) 

where Uj = Uii and Ui is the magnitude of the ith component of the velocity vector in a 
chosen cartesian coordinate system. Eq. when contracted with the projection operator 
5ij/2 and then averaged over k gives the energy spectrum E{k). It is therefore clear that 
the types of correlations contributing to E{k) are those in which i = j and hence i, j 
and the unit wavevector k all lie in the same plane. The corresponding picture in real 
space is to consider the index-symmetric two-point spatial correlation functions -R^(r) = 
|(Mj(x)uj(x -|- r) -|- Uj{^x)ui{x + r)) which has the tensor representation ~ A{r)6ij + B{r)^-^ 
for the isotropic case; the incompressibility constraint gives a relationship between A{r) 
and B{r). This index-symmetric correlation function thus has non-zero contributions when 
i, j and r are co-planar. We will refer to these as 'in-plane' correlations (see Fig. ^ for a 
sketch of these types of isotropic correlations). Similarly, Eq. (0) when contracted with the 
antisymmetric curl operator isijiki, where i = and then averaged over k gives the total 

helicity density H{k) = 2kE{k). (Note that this relationship is distinct from the Schwartz- 
inequality \H{k)\ < 2kE{k)). Therefore, the types of correlations contributing to Eij{k) 
(and hence to H{k)) are those in which i, j and unit wavevector k are mutually orthogonal. 
Again, the corresponding formulation in real space is the index-antisymmetric two-point 
spatial correlation functions -R^(r) = i(Mj(x)nj(x + r) — Mj(x)uj(x + r)) which has the tensor 



representation ~ SijiTi/r and thus has non-zero contributions when i, j and f are mutually 
orthogonal to each other (see Fig. El for a sketch). We will refer to these as 'out-of-plane' 
correlations. Here E{k) = J2\k\=k ^^^^ H{k) = J2\k\=k^0^) ■'*^(~k) where u(k) and 

a;(k) are the fourier transforms of the velocity u(x) and the vorticity uj{'x) = V x u(x) 
respectively. 

The Kraichnan time-scale for energy transfer, te, corresponds to correlations of the type 
Eij{h) (Eq. (P)) which arise due to shearing motions in the plane of coordinates i, j and unit 
wavenumber k jsj . Such in-plane shearing motions cannot give rise to correlations of the type 
£^jj(k) (Eq. (0)) which relate orthogonal components Uji and Uj} across the third mutually 
orthogonal direction k. For this we require out-of-plane shearing motions as depicted in 
Fig. 121 which are provided by the presence of helicity 0,0|- We first derive the time-scale, 
th, associated with such an out-of-plane shear. The governing factor is the relative helicity 
\H{k)\/{2kE{k)) which will be shown to fall off linearly in wavenumber restoring parity as 
k becomes very large. Crucially, we will show that the ratio te/th ~ {\H{k)\/ {2kE{k))Y/^, 
which decays slower than the relative helicity. Therefore, the effect of th cannot be neglected. 
We demonstrate the effect of this new timescale on energy and helicity spectra, and offer 
an interpretation of the "bottleneck" effect observed in measured energy spectra. Finally, 
the new dynamics reveal a dissipation scale which is larger than the Kolmogorov dissipation 
scale, suggesting that the joint cascade is truncated sooner in wavenumber space if helicity 
is present. 

We performed two simulations of the three-dimensional {3-d) , forced Navier-Stokes equa- 
tion in a unit-periodic box with 512 (data I), and 1024 (data II) grid points to a side 
respectively. In these units, the wavenumber k is in integer multiples of 2tt. Energy and 
helicity were iniected into the flow for /c < 2 at each time-step. The forcing scheme was 

n " nil 

the same as in For case I we imposed maximum helicity in < 2 |^ resulting in 
a mean helicity over time of -26.8 in the units of our simulation. For case II the helicity 
input was uncontrolled and random, resulting in a mean helicity of -0.12 which is essentially 
zero compared to case I. The spectra for case I were averaged over 40 snapshots spanning 
8 large-eddy turnover times after spin-up. The spectra for case II were averaged over 48 
snapshots spanning 2 large-eddy turnover times after spin-up. The spin-up time in each case 
was defined to be when the input rate of energy matched the dissipation rate of energy, the 
flow having achieved statistically steady state. Additional parameters of the simulations are 
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given in Table HI 

We recall the introduction in [5] of the distortion timescale (or eddy-turnover time) of an 
eddy with wavenumber k 



k 

Eip)p'dp) ~ [E{k)k') . (3) 

where Kraichnan asssumes that only wavenumbers < k will have a shearing action on 
wavenumbers of order k; the effects from wavenumbers > k will average out. Notice that 
the local shear timescale thus defined depends on the in-plane correlations which contribute 
to the energy spectrum as described above. Analogously, we can define the time-scale th 
for out-of-plane distortions of an eddy, from the antisymmetric co-spectrum, 

^H-{j\m\p'dp)'' ~ {^-\H{k)\kY- (4) 

The distortion or shear corresponding to th is different from the distortion corresponding to 
te of 5j . The shear in the former case derives from the out-of-plane correlations contributing 
to the antisymmetric co-spectrum and hence to the helicity spectrum. The shear motion in 
this case is an out-of-plane twist for which the time-scale is different than the in-plane shear 
rate corresponding to te- 

We estimate the transfer rate (flux) of helicity through wavenumber fc, Fnik) ~ 
k\H{k)\/TH ~ k'^\H(k)f^'^. Assuming steady-state, inertial range behavior with constant 
flux of helicity Fnik) = h, the mean helicity dissipation rate, we obtain 

H{k) ^ h^'^'k-^'K (5) 



1 / \ -1 



We compare with the time-scale for energy transfer of Eq. (jS)), 

TE f\H{k)\\^n 



(6) 



Th \2kE{k)J 

Since [ \H{k)\/{2kE{k)) ]^/^ > \H{k)\/{2kE{k)) as the latter tends to zero, Eq. © implies 
that even for small values of the relative helicity, the time-scales can become comparable. 
This is a fundamental point of difference from previous works in which the presence of 
helicity was considered inconsequential j^. In previous arguments the fact that relative 
helicity must go to zero as 1/ /c in wavenumber space meant that helicity could not have an 
effect on the long term dynamics since it must eventually be dominated by energy, restoring 
parity. Our present analysis shows that while the relative helicity does indeed go rapidly to 



zero, the relative timescale of helicity and energy transfer vanishes much slower. In other 
words, while the energy timescale is always faster, because of the Schwartz equality, the 
helicity timescale can remain comparable to it well into the large wavenumbers. (For the 
initial value problems (decaying turbulence), the evolution equation of the relative helicity 
H{k) / {2kE{k)) and its analytical bounds can be found in 10[.) 

With this second timescale at hand, we can now justifiably ask what the effect of th will 
be on the energy spectrum. The energy flux e through wavenumber k is 

e ^ hE!^ ^ kE{k)H{kY'^k. (7) 
Using Eq. (0) in Eq. ((7j), we get 

E{k) ^ eh-^'''k-^'\ (8) 

To summarize thus far, the te dynamics result in Kolmogorov k~^f^ scaling in both energy 
and helicity spectra, whereas th dynamics result in k~^/^ scaling in both. Clearly, the steeper 
/^-5/3 gcaiing should dominate in the low-wavenumbers while the k~^^^ should manifest in 
the higher wavenumbers. We emphasize that in order for the latter scaling to be visible in 
the high wavenumbers, th cannot be too much slower than te- As shown above, according to 
Eq. (jH)), this may occur for very modest relative helicity in the high wavenumbers, contrary to 
previous assumptions. The main point of our paper is that in general, the helicity timescale 
Th niay not be ignored. 

Figure El shows the energy and helicity spectra from our simulations. Each of the spectra 
are compensated by fc^^^ and by k^^^ in order to distinguish the dominant scaling. For the 
strongly helical case I, there is good agreement with k~^/^ scaling for both the energy and 
helicity spectra in approximately the same range (Figures|2^a) and (b)) of slightly less than a 
decade. Note that the compensation with k^/^ results in the commonly observed 'bottleneck' 
phenomenon which we will discuss below. The energy spectrum of II (Fig. Efc)) shows a 
range of k~^^^ scaling followed by a range of k~^^^ scaling (which appears as a bottleneck 
in the k~^/^ compensated plot). The scaling ranges are modest even at this high resolution 
of 1024^, but nonetheless the results are telling; the scaling is most certainly not k~^/^ 
throughout and the agreement with k~^^^ though over a short range, certainly indicates 
shallower than k~^/'^ scaling behavior of the bottleneck region. The relative helicities in 
the range 10 < /c < 100, where the k~^^'^ scaling is seen, are shown in Fig. |3] For I, the 



relative helicity falls from about 10% to about 3% corresponding to te/th ranging from 32% 
to 17% according to Eq. ©• Despite the negligibly small total helicity H = —0.12 of II, 
and its noisy helicity spectrum (not shown), its relative helicity values lie between 1% and 
5%. This implies that te/th could be as much as 22%. In both cases th might in fact not 
be much longer than te- It is important at this stage to comment on the appearance of a 
helicity-dependent scaling feature in flow II which is nominally helicity-free on average. The 
flrst point is that zero average helicity does not imply that the average helicity spectrum 
H{k) is zero for all k. In fact, we only input energy and helicity in the low-wavenumbers, 
the Navier-Stokes dynamics determines the helicity in all other wavenumber, including the 
highest wavenumbers where in fact there is the well-known viscous helicity production. 
There is therefore no control of the helicity in an given wavenumber and the spectrum is 
generally not zero everywhere. The second point is that our analysis shows that it is not 
the total helicity but the relative helicity which determines the trade-off between the two 
timescales. Given these two points it is not contradictory to measure k~^^^ spectral scaling 
in the flow with negligible mean helicity. In fact, this flow is probably more similar to most 
experimental flows which are close to helicity-free in the mean but with uncontrolled and 
often unknown helicity spectra 11 1. 

These results are the flrst indication of the possibility of k~'^^^ scaling ranges in both 
energy and helicity spectra simultaneously. The possibility of a 'pure' or 'maximal' forward 
cascade of helicity scaling as k'"^^^ Jd], with inverse cascade of energy scaling as does 
not arise. This is because in our analysis we have retained the effect of the helical time-scale 
th and allowed it to modify the spectral dynamics. Since the scaling corresponding to th is 
k~^^^, a slower decay than k~^^'^, its 'signature' in the spectra can dominate at large k even 
as the overall parity is being restored. 

Based on the analysis above, we propose that the bottleneck in the total energy spectrum 
is in fact a change in the scaling of the energy spectrum, from a k~^^^ regime in which the te 
dynamics dominate to a less steep k~^^^ regime in which the th dynamics become significant. 
We will use the kinematic arguments of [1^ with our new phenomenology and dynamics to 
analyse the bottleneck for such a helical influence. In simulations, it is possible to compute 
the total energy spectrum E{k) = (1/2) X]|k|=fe ^ '^^^^ ^ ^^^^^ radius k rather 

accurately. In experiments it is convenient to measure the one-dimensional (l-rf) longitudinal 
and transverse spectra along the measurement direction, say z. In our 3-d flow simulation 
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we calculate the 1-d spectra as follows. The 1-d fourier transform in the 2;-direction of 
the velocity u(x) is u{x,y,kz) = (1/A^) ^^^^^ e*'^^^"u(x, z„) where < A;^ < tt/6z. The 
longitudinal (transverse) 1-d spectrum, averaged over the x-y plane, is defined by 

1 ^ 

EL{T){kz) = j^j^ Yl \^zi±}{xp,yg,k,)\'^, (9) 

where \u±\'^ = \ux\'^ + In isotropic flow, the 1-d spectra should be independent of 

the direction in which the fourier transform is performed. Our time-averaged longitudinal 
spectra computed in the the x, y and z directions all collapse and the transverse spectra do 
the same. We average the 1-d spectra in the three coordinate directions and drop the use 
of the subscript z to denote the direction of the fourier transform. In isotropic flow there is 



a relation between the 1-d and 3-d spectra 
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3, 



Our spectra satisfy Eq. ()10|) very well for A; > 10 and fairly well in the lower wavenumbers. 
As emphasized in 1^, Eq. (fTUI) is a local relationship, wherein the functional form of the 
total spectrum E{}z) is fully determined by the local behavior of Ei^ and Ej^ at a given wave 
number. For homogeneous, isotropic flows with helicity, it is reasonable to suppose that 
Eilk) (see definition in Eq. (jHl)) would mainly carry contributions from the in-plane shear 
time-scale te- However Ex{k), which is related to the transverse components of the velocity 
fourier transform, could be influenced by th dynamics coming from E{k). The correlation 
time between transverse components can be slowed down by the dynamics of E{k) which 
arise due to the presence of helicity. Such coupling may not be deduced from kinematic argu- 
ments, it requires proper consideration of the new dynamics. Furthermore, it is not possible 
to see this coupling in the unclosed lowest-order Karman-Howarth dyanamical equations 
wherein the symmetric and antisymmetric parts completely decouple for homogeneous flows 
I?]. It was pointed out in that higher-order statistical equations for helical turbulence 
should show the coupling. The coupling is seen explicitly in the EDQNM closure [l^ . 

Figure El shows the compensated longitudinal and transverse spectra for simulation II. 
The bottleneck is greatly diminished in Ei{k) (Fig. El^a), solid curve) which shows close to 
^-5/3 gcaiing throughout the inertial range. In Exik) (Fig. Efb), solid curve), the bottle- 
neck persists although its peak occurs slightly earlier in wavenumber than the bottleneck 
for the corresponding total spectrum. For completeness, we have also shown the 1-d spectra 



compensated by /c^/^ (Fig. El dotted curves). Based on our arguments above, we might 
have expected a stronger k~^^^ scahng region of the transverse spectrum; such behavior is 
not clearly observed although there is a tendency towards a scaling shallower than k~^^^ in 
the bottleneck regime. We plan in future work to check the present indications that the 
bottleneck will be stronger in the transverse spectrum than in the longitudinal one because 
of the greater contamination of the former by the helical co-spectrum dynamics. This raises 
the intriguing possibility that th affects the scaling of transverse structure functions, ac- 
counting for some of the observed difference between the scaling exponents of longitudinal 
and transverse structure functions in near-isotropic, high-Reynolds number turbulence data 
[isj ]. Such contributions would appear as a parity-breaking in the isotropic small-scales, 
and might not be easily disentangled by, for example, the SO (3) group decomposition meth- 
ods (see lis] and references therein) used to extract isotropic contributions to non-helical 
turbulence statistics. Although this aspect of the influence of helicity dynamics remains 
speculative we hope this work provides sufficient motivation for further investigation. 

The bottleneck is a well-known phenomenon which has been observed in experimental 
measurements 17,^^ and Navier-Stokes simulations jiol. I20I1. While the mean helicity of 
the flows in these investigations might be zero (although many do not report the mean 
helicity, using the often reasonable assumption that the flow is non-helical on average), their 
relative helicities might not be, and indeed have not been reported, because the connection 
between the Kolmogorov phenomenology, helicity dynamics and the bottleneck did not exist. 
There have been various different approaches taken to explain this phenomenon, including 



viscous effects 



2l| and various kinematic arguments used to fit to a parametrized form 



22L I23L |2J|. To the best of our knowledge, there has thus far never been an association of 
the bottleneck with helicity. What we are proposing here is a fundamental physical cause of 
the bottleneck due to the helicity dynamics slowing down the cascade of energy and helicity, 
the two conserved quantities in turbulence. Further, our empirical evidence, particularly 
in case II, where the total helicity is negligibly small, indicates that this effect could occur 
even in flows with essentially zero mean helicity but with non-zero relative helicity spectra 
(i.e. H{k) is not zero everywhere). Said differently, even if j H{k)dk = 0, there can be a 
range of k where \H{k)\ and hence \H{k)\/2kE{k) is finite and possibly large enough that 
te/th ~ {\H{k)\/2kE{k)y/^ is not negligible. This effect was reported in an experimental 
work investigating spontaneous reflection-symmetry breaking in boundary layer flows 
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Such a scenario may occur since, while global helicity is statistically conserved in the inertial 
range, local helicity is not. Let ^(x) be the local helicity ^(x) = ■u(x) ■ (^{x). The equation 
for ^ reads 

1 



dt^ + Ujdji = -djiujjp) + -dj{ujj\u 



2\ 



+v[V^^-2{diU,){diu:j)]. (11) 

Locally, both the nonlinear and the viscous terms in the helicity dynamics might play a role 
in enhancing or diminishing the helicity. If the non-zero relative helicity in the scaling range 
arises from the nonlinear term, then the effects we see are indeed valid in the high- Reynolds 
number (inviscid) limit; if they arise from the viscous term then the effects we see would 
disappear at very high Reynolds numbers. This hypothesis is not testable at the present 



time but our work, which studies data from simu 
several performed before (see for example ^, 25 



ations which are similar to, or the same as. 



26|), with comparable Reynolds number. 



gives substantial motivation to further examine these questions. We hope in particular 
to motivate measurements of the relative helicity and scaling behavior of the bottleneck 
region, in other flows which report the bottleneck phenomenon in order to further check this 
connection. 

We finally present a key result from analysis of the convergence of the dissipation integrals 



for a two-timescale cascade. In 



Q 



the same analysis was performed assuming a single 
timescale te- The total dissipation of energy and helicity may be written as 

De = 2u{^ f ^ dkP E{k) + f " dkk"^ E{k)^ (12) 

Dh = 2z/(^^ ^ dk k^ H{k) + ' dk k^ H{k)^, (13) 

where 

E{k) = <^ (14) 
I Ch eh-^/'^k-^/^ for k^ < k < kd, 

{CEE^^/^k"^/^" ioik<kc\ 
(15) 
Ch /i2/^A;-4/3 for k^ < k < kd, 

and Cei Ch-, ce and ch are constants. The precise estimate of the transition wavenumber kc 
is unimportant to what follows. The upper wavenumber kd denotes the maximum beyond 
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which the cascade is completely suppressed by viscosity. We obtain 

De-^u e h-^/'^ k]/^ (16) 
Dur^v h^/'kf. (17) 

In ^ kd = ke (£/z/3)V4 (the Kolmo gorov dissipation wavenumber). In our case, set- 
ting kd = ke causes the integrals to diverge as z/^^/^ in the limit ^ 0. The choice 
kd ~ {h/v'^Y^^ ensures that the integrals converge to their correct values for statistically 
steady-state, De = £ and Dh = h. We will call this new wavenumber kh since it depends 
solely on the helicity dissipation rate. It must be distinguished from the kn of which 
depends on both energy and helicity dissipation rates. In the limit u —>■ 0, k^ ^ kh. In our 
simulations k^ > k^hj a factor of about 2.5. While in agreement with our analysis, we can- 
not really distinguish between the two wavenumbers in this data. However, we suggest that 
the resolution requirement for measurements in turbulence with helicity, or more precisely, 
with non-zero helicity spectra, might be weaker than that in turbulence without helicity. 
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FIG. 1: The in-plane longitudinal and transverse correlation configurations which contribute to 
the isotropic symmetric correlation function Rfj{r). 



u,(x) 




FIG. 2: The out-of-plane correlation configuration which contributes to the isotropic antisymmetric 
correlation function R^{r). The intrinsic 'handed 
arrow, cannot appear in the geometry of Fig. ^ 



correlation function R^{r). The intrinsic 'handedness' of this configuration, indicated by the curved 
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FIG. 3: Compensated spectra of (a) energy and (b) helicity, for data I. Observe the k~^^^ scaling 
(dashed curve) range indicated by the horizontal line segment, (c) Compensated energy spectrum 
for data II. The horizontal line segments indicate ranges of k^^^^ (solid curve) and k~^^^ (dashed 
curve) scaling. 
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FIG. 4: The relative helicities for (a) data I and (b) data II, in the range 10 < A; < 100, where 
Fig. 13] shows k~'^^^ scahng. The relative helicity lies between 3% and 10% in (a) and between 1% 
and 5% in (b). The corresponding ratio of timescales te/th (Eq. @), is estimated to be upto 32% 
for data I upto 22% for data II. 
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FIG. 5: Compensated (a) longitudinal and (b) transverse energy spectra for the simulation II. 
Ei{k) shows a diminished bottleneck where the full spectrum E[k) of Fig. Ofc) has a pronounced 
bottleneck. Exik) shows a range of k~^/'^ coinciding with that of Figl^fc) and a persistent bottle- 
neck with a scaling shallower than k~^^^ scaling in the region of the bottleneck. 



Tables 





N 


u X 10^ 


Rx 


E 


e 


H 


h 


r]e X 10'^ 


Vh X 10^ 


I 


512 


1 


270 


1.72 


1.51 


-26.8 


62.2 


1.7 
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II 


1024 


0.35 


430 


1.87 


1.75 


-0.12 


13.2 


1.3 


4 



TABLE I: Parameters of the numerical simulations I and II. u - viscosity; Rx - Taylor Reynolds 
number; mean total energy E = ^ Ylk ^ ' mean energy dissipation rate; mean total helicity 

H = Y^j^ u(k) • ci;(— k); h - mean helicity dissipation rate; = (u^/e)^/^; r]h = {u^/h)^/^. 
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